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Abstract. In this communication, complex systems with a near trivial 
dynamics are addressed. First, under the hypothesis of equiprobability 
in the asymptotic equilibrium, it is shown that the (hyper) planar geom- 
etry of an iV-dimensional multi-agent economic system implies the ex- 
ponential (Boltzmann-Gibss) wealth distribution and that the spherical 
geometry of a gas of particles implies the Gaussian (Maxwellian) dis- 
tribution of velocities. Moreover, two non-linear models are proposed to 
explain the decay of these statistical systems from an out-of-equilibrium 
situation toward their asymptotic equilibrium states. 
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1 Introduction 

In this paper, different classical results [112] are recalled. They are obtained from 
a geometrical interpretation of different multi-agent systems evolving in phase 
space under the hypothesis of equiprobability [314] . Two nonlinear models that 
explain the decay of these statistical systems to their asymptotic equilibrium 
states are also collected [516] . 

We sketch in section [5] the derivation of the Boltzmann-Gibbs (exponen- 
tial) distribution [3] by means of the geometrical properties of the volume of 
an iV-dimensional pyramid. The same result is obtained when the calculation 
is performed over the surface of a such ./V-dimensional body. In both cases, the 
motivation is a multi-agent economic system with an open or closed economy, 
respectively. 

Also, a continuous version of an homogeneous economic gas-like model [5] is 
given in the section[5] This model explains the appearance, independently of the 
initial wealth distribution given to the system, of the exponential (Boltzmann- 
Gibbs) distribution as the asymptotic equilibrium in random markets, and in 
general in many other natural phenomena with the same type of interactions. 

The Maxwellian (Gaussian) distribution is derived in section [3] from geo- 
metrical arguments over the volume or the surface of an iV-sphere [3]. Here, 
the motivation is a multi-particle gas system in contact with a heat reservoir 
(non-isolated or open system) or with a fixed energy (isolated or closed system), 
respectively. 



The ubiquity of the Maxwellian velocity distribution in ideal gases is also 
explained in the section[3]with a nonlinear mapping acting in the space of velocity 
distributions [B] . This mapping is an operator that gives account of the decay of 
any initial velocity distribution toward the Gaussian (Maxwellian) distribution. 

Last section contains the conclusions. 



2 Systems Showing the Boltzmann-Gibbs Distribution 
2.1 Multi- Agent Economic Open Systems 

Here we assume N agents, each one with coordinate Xj, i = 1, . . . ,N, with x, > 
representing the wealth or money of the agent i, and a total available amount 
of money E: 

xi + x 2 + • • • + xjv-i + xn < E. (1) 

Under random or deterministic evolution rules for the exchanging of money 
among agents, let us suppose that this system evolves in the interior of the 
./V-dimensional pyramid given by Eq. (fTJ). The role of a heat reservoir, that in 
this model supplies money instead of energy, could be played by the state or 
by the bank system in western societies. The formula for the volume Vn(E) of 
an equilateral ./V-dimensional pyramid formed by N + 1 vertices linked by N 
perpendicular sides of length E is 

Vn{E) = — . (2) 

We suppose that each point on the ./V-dimensional pyramid is equiprobable, then 
the probability / (x,)dxj of finding the agent i with money Xi is proportional to 
the volume formed by all the points into the [N— l)-dimensional pyramid having 
the ith-coordinate equal to Xj. Then, it can be shown that the Boltzmann factor 
(or the Maxwell-Boltzmann distribution), /(Xj), is given by 

that verifies the normalization condition 

fE 

f(xi)dxi = 1. (4) 







The final form of /(x), in the asymptotic regime N — > oo (which implies E — > oo) 
and taking the mean wealth e = E/N, is: 

/(x)rfx = - e~ x / e dx, (5) 
e 

where the index i has been removed because the distribution is the same for 
each agent, and thus the wealth distribution can be obtained by averaging over 
all the agents. This distribution has been found to fit the real distribution of 
incomes in western societies [7]. 



2.2 Multi- Agent Economic Closed Systems 



We derive now the Boltzmann-Gibbs distribution by considering the system in 
isolation, that is, a closed economy. Without loss of generality, let us assume 
N interacting economic agents, each one with coordinate Xi, i = 1, . . . , N, with 
Xi > 0, and where Xi represents an amount of money. If we suppose that the 
total amount of money E is conserved, 

x\ + x 2 + • • • + Xn-x + x N — E, (6) 

then this isolated system evolves on the positive part of an equilateral N- 
hyperplane. The surface area Sn(E) of an equilateral ./V-hyperplane of side E is 
given by 

S N (E) = ^ E»-K (7) 

If the crgodic hypothesis is assumed, each point on the iV-hyperplane is equiprob- 
able. Then the probability f(xi)dxi of finding agent i with money Xi is propor- 
tional to the surface area formed by all the points on the iV-hyperplane having 
the ith-coordinate equal to Xi. It can be shown that f{x{) is the Boltzmann 
factor (Boltzmann-Gibbs distribution), with the normalization condition (|4|. It 
takes the form, 

ft \ 1 S N -x{E-Xi) 

f{Xi) = — , (8 

Sn(E) smffN 



where the coordinate 9n satisfies sin^jv = y^jr~- After some calculation the 
Boltzmann distribution is newly recovered: 

f{x)dx = — e- x/kT dx, (9) 

KT 

with e = kr, being k the Boltzmann constant and r the temperature of the 
statistical system. 



2.3 The Continuous Economic Gas-like Model 

We consider an ensemble of economic agents trading with money in a random 
manner [7]. This is one of the simplest gas- like models, in which an initial amount 
of money is given to each agent, let us suppose the same to each one. Then, pairs 
of agents are randomly chosen and they exchange their money also in a random 
way. When the gas evolves under these conditions, the exponential distribution 
appears as the asymptotic wealth distribution. In this model, the microdynamics 
is conservative because the local interactions conserve the money. Hence, the 
macrodynamics is also conservative and the total amount of money is constant 
in time. 

The discrete version of this model is as follows [7] . The trading rules for each 
interacting pair (m^ rrij) of the ensemble of N economic agents can be written 



as 



m- = a (mi + rrij), 

rrij = (1 — <j)(rrii + rrij), (10) 

i-.-.v. 

where a is a random number in the interval (0,1). The agents (i, j) are randomly 
chosen. Their initial money (mi, rrij), at time t, is transformed after the interac- 
tion in (m^m'j) at time t + 1. The asymptotic distribution pf(m), obtained by 
numerical simulations, is the exponential (Boltzmann-Gibbs) distribution, 

Pf(m) — /3exp(-/3m), with /3 = 1/ < m > gas , (11) 

where pf(m)dm denotes the PDF (probability density function), i.e. the proba- 
bility of finding an agent with money (or energy in a gas system) between m 
and m + dm. Evidently, this PDF is normalized, = J pf(m)dm = 1. The 
mean value of the wealth, < m > ga s, can be easily calculated directly from the 
gas by < m > gas = J2i mjN. 

The continuous version of this model [5] considers the evolution of an initial 
wealth distribution po(m) at each time step n under the action of an operator 
T . Thus, the system evolves from time n to time n + 1 to asymptotically reach 
the equilibrium distribution pf(m), i.e. 

lim T n (p (m))^p f (m). (12) 

n— >oo 

In this particular case, pf(m) is the exponential distribution with the same 
average value < pf > than the initial one < po >, due to the local and total 
richness conservation. 

The derivation of the operator T is as follows [5]. Suppose that p n is the 
wealth distribution in the ensemble at time n. The probability to have a quantity 
of money x at time n + 1 will be the sum of the probabilities of all those pairs 
of agents (u, v) able to produce the quantity x after their interaction, that is, 
all the pairs verifying u + v > x. Thus, the probability that two of these agents 
with money (u,v) interact between them is p n (u) * p n (v). Their exchange is 
totally random and then they can give rise with equal probability to any value 
x comprised in the interval (0,it + v). Therefore, the probability to obtain a 
particular x (with x < u + v) for the interacting pair (u,v) will be p n (u) * 
p n (v)/(u + v). Then, T has the form of a nonlinear integral operator, 

t \ T f \ f f Pn(u)Pn(v) , , 

p„+i(x) = Tp n (x) = // dudv. (13) 

J Ju+v>x U + V 

If we suppose T acting in the PDFs space, it has been proved [5] that T 
conserves the mean wealth of the system, < Tp >=< p >. It also conserves 
the norm (|| • ||), i.e. T maintains the total number of agents of the system, 
\\Tp\ \ = \ \p\ \ — 1, that by extension implies the conservation of the total richness 
of the system. We have also shown that the exponential distribution p / (x) with 



the right average value is the only steady state of T, i.e. Tpf = pf. Computations 
also seem to suggest that other high period orbits do not exist. In consequence, 
it can be argued that the relation (IT^I) is true. This decaying behavior toward the 
exponential distribution is essentially maintained in the extension of this model 
for more general random markets. 



3 Systems Showing the Maxwellian Distribution 
3.1 Multi-particle open systems 

Let us suppose a one-dimensional ideal gas of N non-identical classical particles 
with masses mi, with i = 1, . . . , N, and total maximum energy E. If particle i 
has a momentum rriiVi, we define a kinetic energy: 

Ki = p 2 i = ^mivl (14) 

where pi is the square root of the kinetic energy JQ. If the total maximum energy 
is defined as E = R 2 , we have 

v\+pI + ---+Pn-i+P 2 n<R 2 - (15) 

We see that the system has accessible states with different energy, which can be 
supplied by a heat reservoir. These states are all those enclosed into the volume 
of the iV-sphere given by Eq. (|15l) . The formula for the volume Vn{R) of an 
iV-sphere of radius R is 

JV 

Vlf W = r(l*+i) RN ' (16) 

where -T(-) is the gamma function. If we suppose that each point into the iV- 
sphere is equiprobable, then the probability f(pi)dpi of finding the particle i 
with coordinate pi (energy pf) is proportional to the volume formed by all the 
points on the iV-sphere having the ith-coordinatc equal to pi. It can be shown 
that 

'<*> - Y ^Sv Sl - 

which is normalized, f_ R f(Pi)dpi = 1. The Maxwellian distribution is obtained 
in the asymptotic regime N — > oo (which implies E — > oo): 



f(p)dp=^ e e-^^dp, (18) 

with e — E/N being the mean energy per particle and where the index i has 
been removed because the distribution is the same for each particle. Then the 
equilibrium velocity distribution can also be obtained by averaging over all the 
particles. 



3.2 Multi-particle closed systems 



We start by assuming a one-dimensional ideal gas of N non-identical classical 
particles with masses m^, with i = 1, . . . ,N, and total energy E. If particle i 
has a momentum rriiVi, newly we define a kinetic energy Ki given by Eq. (1141) . 
where pi is the square root of Ki. If the total energy is defined as E = R 2 , we 
have 

pI+P 2 2 + ---+P 2 n-i+P 2 n = R 2 - (19) 

We see that the isolated system evolves on the surface of an TV-sphere. The 
formula for the surface area Sn(R) of an A^-sphere of radius R is 



27rt 

r(f)' 



S N (R) = T^lvT^" 1 ' (2°) 



where r(-) is the gamma function. If the ergodic hypothesis is assumed, that 
is, each point on the ./V-sphere is equiprobable, then the probability f(pi)dpi of 
finding the particle i with coordinate Pi (energy p 2 ) is proportional to the surface 
area formed by all the points on the iV-sphere having the ith-coordinate equal 
to pi . It can be shown that 



f(p) - 1 gwVgEl) (21) 

m) ~S N (R) (1-11)1/2 ' (21) 

which is normalized. Replacing p 2 by ^mv 2 , f(p) takes the following form g(v) 
in the asymptotic limit N — > oo, 



9{v)dv = J^-e-™*l 2kT dv. (22) 
V 2ttkt 

This is the typical form of the Maxwellian distribution, with e = fer/2 given by 
the equipartition theorem. 



3.3 The Continuous Model for Ideal Gases 

Here, as we have done in the anterior case of economic systems, we present a 
new model to explain the Maxwellian distribution as a limit point in the space 
of velocity distributions for a gas system evolving from any initial condition [B] . 

Consider an ideal gas with particles of unity mass in the three-dimensional 
(31?) space. As long as there is not a privileged direction in the equilibrium, 
we can take any direction in the space and study the discrete time evolution of 
the velocity distribution in that direction. Let us call this direction U. We can 
complete a Cartesian system with two additional orthogonal directions V, W. If 
p n (u)du represents the probability of finding a particle of the gas with velocity 
component in the direction U comprised between u and u + du at time n, then 
the probability to have at this time n a particle with a 3D velocity (u, v, w) will 
be p n (u)p n {v)p n (w). 



The particles of the gas collide between them, and after a number of interac- 
tions of the order of system size, a new velocity distribution is attained at time 
n + 1. Concerning the interaction of particles with the bulk of the gas, we make 
two simplistic and realistic assumptions in order to obtain the probability of hav- 
ing a velocity x in the direction U at time n+ 1: (1) Only those particles with an 
energy bigger than x 2 at time n can contribute to this velocity x in the direction 
U, that is, all those particles whose velocities (u, v, w) verify u 2 + v 2 + w 2 > x 2 ; 
(2) The new velocities after collisions are equally distributed in their permitted 
ranges, that is, particles with velocity (u,v,w) can generate maximal velocities 
±Umax = ±V" 2 + v 2 + w 2 , then the allowed range of velocities [—U max , U max ] 
measures 2|{/ ma;E |, and the contributing probability of these particles to the ve- 
locity x will be Pn{u)pn(v)p n (w) / {2\U m ax\)- Taking all together we finally get 
the expression for the evolution operator T. This is: 

Pn{u)p„(v)p n (w) 



Pn+i{x) = Tp n {x) = / / / ==■ dudvdw . (23) 

lu 2 +v 2 +w 2 >w 2 2yu + v z + w z 



Let us remark that we have not made any supposition about the type of 
interactions or collisions between the particles and, in some way, the equivalent 
of the Boltzmann hypothesis of molecular chaos [5] would be the two simplistic 
assumptions we have stated on the interaction of particles with the bulk of 
the gas. Then, an alternative framework than those usually presented in the 
literature |10] appears now on the scene. In fact, it is possible to show that 
the operator T conserves in time the energy and the null momentum of the 
gas. Moreover, for any initial velocity distribution, the system tends towards its 
equilibrium, i.e. towards the Maxwellian velocity distribution (ID case). This 
means that 

lim T" (po(x)) -> Pa (x) = \[^e~ ax2 (24) 

n— >oc y 7T 

with a — (2 < x 2 ,po >) . In physical terms, it means that for any initial 
velocity distribution of the gas, it decays to the Maxwellian distribution, which 
is just the fixed point of the dynamics. Recalling that in the equilibrium < 
x 2 ,p a >= fcr, with k the Boltzmann constant and r the temperature of the gas, 
and introducing the mass m of the particles, let us observe that the Maxwellian 
velocity distribution can be recovered in its 3D format: 

p a {u) Pa {v)p a {w) = (— )* e -™*(« 2 +«W) with a = (2k T )- 1 . (25) 

In general, it is observed that the convergence of the T-iterations of any distri- 
bution p(x) to its Gaussian limit p a (x) is very fast. 



4 Conclusion 



We have shown that the Boltzmann factor describes the general statistical be- 
havior of each small part of a multi-component system whose components or 



parts are given by a set of random variables that satisfy an additive constraint, 
in the form of a conservation law (closed systems) or in the form of an upper 
limit (open systems). 

Let us remark that these calculations do not need the knowledge of the exact 
or microscopic randomization mechanisms of the multi-agent system in order to 
reach the equiprobability. In some cases, it can be reached by random forces, in 
other cases by chaotic or deterministic causes. Evidently, the proof that these 
mechanisms generate equiprobability is not a trivial task and it remains as a 
typical challenge in this kind of problems. 

In order to explain the ubiquity and stability of this type of distributions 
two models based on discrete mappings in the space of distributions have been 
proposed. On one hand, the gas-like models interpret economic exchanges of 
money between agents similarly to collisions in a gas where particles share their 
energy. The continuous version of a gas-like discrete model where the agents 
trade in binary collisions has been introduced to explain the stability of the 
exponential distribution in this kind of economic systems. On the other hand, 
a nonlinear map acting on the velocity distribution space of ideal gases, which 
gives account of the decay of an out-of-equilibrium velocity distribution toward 
the Maxwellian distribution, has been presented. Some properties concerning the 
dynamical behavior of both operators have also been sketched. 
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